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1 Introduction 

The integration of systems of differential equations admitting infinitesimal symmetries 
was the main concern of Lie in developing what is nowadays called the theory of Lie 
algebras and Lie groups. In particular, in a remarkable work ^5] he was able to prove 
an important theorem connecting Lie algebras and nonlinear superposition rules for 
solutions of some non-autonomous systems of nonlinear ordinary differential equations. 
These systems can be considered as a generalization of linear systems but the superpo- 
sition rule is no longer a linear function. Our aim in this paper is to study once again, 
from a geometric viewpoint, the theory of systems of differential equations admitting 
a (maybe nonlinear) superposition rule, allowing us to express the general solution of 
the system by means of a superposition function in terms of a (fundamental) set of 
particular solutions, with the hope of establishing clearly the necessary and sufficient 
conditions for a system to admit such a superposition rule. This time, however, we 
include also partial differential equations into our considerations. 

Even if the hypotheses of Lie's theorem were not accurately stated from the today 
level of rigor, the resulting systems characterized by means of an associated Lie algebra 
appear very often in physical problems and in many cases the problem is related with 
another one on the corresponding Lie group. This provides us with both methods of 
reduction to simpler problems on one side, and another method, introduced by Wei 
and Norman which involves some algebraic manipulations based on Lie groups and Lie 
algebra theories, on the other. 

As far as we know, there is no rigorous proof of the if part in Lie's theorem and 
the attempts known to us to get a rigorous geometric proof share the same pseudo- 
argument [21 E] • In this paper we prove that actually the existence of a superposition 
rule for the solutions of a given non-autonomous system implies that it has the explicit 
form which is usually accepted. The proof is based on an alternative but equivalent 
definition of superposition rule: we consider it as a foliation with some appropriate 
properties explicitly formulated later on. An auxiliary lemma is necessary to overcome 
the weak point in previous derivations of the theorem |2~j I17j. 

On the other hand, the converse part is not given in its full generality and almost 
nothing is said about the uniqueness of the superposition function for these Lie systems. 
Only in ^H] an example given in |15| . for which there are two different superposition 
functions, is pointed out. Our approach provides us with an answer to this important 
question, the key point being the codimension of the foliation constructed from the 
given Lie algebra of vector fields. Moreover, this codimension is very relevant when the 
action of the Lie algebra of vector fields on the initial manifold is not transitive. 

Finally, as a more convincing argument supporting the use of this alternative defi- 
nition of superposition rule, it will be shown that this definition allows an immediate 
generalization of Lie's Theorem for the case of systems of partial differential equations. 
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The organization of the paper is as follows. Next section discuss the concept of 
superposition function and gives a geometric characterization of such superposition in 
terms of a foliation. Section 3 is devoted to a complete proof of the statement of 
Lie's theorem by establishing a lemma which allows us to overcome the weak point of 
other previous derivations. The number of solutions in a fundamental set is discussed 
in Section 4 and in Section 5 the problem of uniqueness of the superposition rule is 
studied. Lie systems on Lie groups and homogeneous spaces are considered in Section 6 
as most important examples. Finally, a generalization of the Lie's Theorem for the case 
of systems of first-order partial differential equations is given in Section 7. An outlook 
with future applications is given in the last section of the paper. 



2 Superposition rules for ordinary differential equa- 
tions 

By a superposition rule (or a superposition principle) for a given system of ordinary 
differential equations 

fix 1 

— = Y*(t,x), i = l,...,n, (1) 

one usually understands, after [T^], a superposition function $ : M n ( m+1 ) — > W 1 given 
by 

x = $0(i), • • • ,X( m ); ki, . . . , k n ) , (2) 
such that the general solution can be written, at least for sufficiently small t, as 

x(t) =®(x (1 )(t),...,x {m) (t);ki,...,k n ) , (3) 

with {x( a )(t) | a = l,...,m} being a fundamental set of particular solutions of the 
system JT} and k = (ki, . . . , k n ) being a set of n arbitrary constants associated with 
each particular solution. 

The standard example is the system of linear differential equations 

— =5>%.(t)^, i = l,...,n, (4) 
which admits the superposition function, with m = n, 

n 

x = $(x ( i),...,X( n) ;/ci,...,/c n ) = ^ kj X(i) . 

t=i 

Of course, we can obtain every solution by superposing xn)(t), . . . ,xi m ){t), for certain 
ki, . . . ,k n , only if the functions xm(t), . . . , xr m ^(t), are appropriately independent, i.e. 
if they form a fundamental set of solutions. In the above example it means that if the 
matrix 

X(t) = (x\ j) (t))} i,j = l,...,n, (5) 

is invert ible for small t. 
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The order in which the particular solutions are chosen is irrelevant and therefore 
the superposition function should be such that a permutation of two arguments only 
amounts to a change of the parameters k. Note also that it is assumed that the super- 
position function $ does not depend explicitly on the independent variable t, and this 
fact has strong consequences (see later the Lie theorem). 

From a geometric perspective, systems of differential equations as (J1J) appear as 
those determining the integral curves of a t-dependent vector field in IR n , 

i=l 

the generalization to the case of a n-dimensional manifold being immediate. Note that 
then in any point x G N, the t-dependent vector field Y in N, determines not only one 
vector but a linear subspace, spanned by the set of vectors {Y(t,x) \ t G M}, of the 
corresponding tangent space. Actually, under a time re-parametrization the vectors are 
rescaled and, when changing the value of t, different vectors are obtained. In this way 
it defines a 'generalized' distribution for which the dimension of the linear subspace can 
change from one point to another. We will see later on that the case we are interested 
in is such that the distribution defined by the t-dependent vector field is involutive. 

In order to look for superposition rules we need a more geometric picture. Let us 
first observe that, as a consequence of the Implicit Function Theorem, the function 
$(x(i), . . . , X( m ); ■ ) : M. n — > MJ 1 can be, at least locally around generic points, inverted, 
so we can write 

k = ^(£(0), • • .,a?(m)) (6) 
for a certain function \1/ : R n ( m+1 ) — » W 1 . Hereafter in order to handle a short notation 
we start writing x(p\ instead of x. The foliation defined by the function \1/ is now 
invariant under permutations of the (m + 1) variables. 

With some abuse of terminology we will also call the function \1/ a superposition 
function. The relation between <3> and \l/ is given by: 

k = . . . , X( TO )J k U . . . , kn),X(i), . . . , X( m) ) . (7) 

For instance, for the system (J3J) we have 

* (z(o)(f), • • -,X( n )(t)) = X^^x^it), 

where X(t) is the matrix given in (jSJ). This example indicates the obvious fact that, 
in general, the superposition function \l/ is defined on an open dense subset of IR n(m+1 ) 
rather than the whole M n ( m+1 ). 

The fundamental property of the superposition function \1/ is that as 

k = ^(X( )(t), £(!)(£), . . . ,X( m )(t)) , (8) 

the function \l/(x(o), . . . , X( m )) is constant on any (m + l)-tuple of solutions of the system 
(|H). This property is true for any choice of (m + 1) solutions and this means that the 
foliation is invariant under the permutation of the (m + 1) arguments of the function 
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After differentiation of relation (JHJ) with respect to t, as the functions Xi(t) are 
solutions of (HJ, we get 

DV(Y(t,x {0) ),...,Y(t,x (m) )) = 0, (9) 

i.e. 

i=l a=0 ( a ) 

and therefore the 'diagonal prolongations' Y(t, xto), • • • , X( m )) °f t ne ^-dependent vector 
field Y(t,x), given by 

m 

Y(t, x (0 ), • • • , x (m) ) = ^ ' t G K , 

a=0 

where 

n 

y a (t,x (a) ) = ^r(t,x (a) )-— (io) 

i=l (a) 

are t-dependent vector fields on R n x ■ • • x R™ ((m + 1) factors) which are tangent to 
the level sets of \l/ as displayed by Q . 

The level sets of \& corresponding to regular values define a n-codimensional foliation 
T on an open dense subset U C lR n x • • • x W 1 ((m + 1) factors) and the family 
{Y(t), t G M} of vector fields in R n ( m+1 ) consists of vector fields tangent to the leaves 
of this foliation. 

This foliation has another important property. Since on the level set correspond- 
ing to k — (ki, . . . , k n ) G R n and given (xm, . . . , sc( m )) G M nm , there is a unique point 
(>(o) ) G .F fc , namely, ($(x ( i), . . . , x ( 

m) ! ™ J ) (1) j • • • j (to) ) g ^ fc (cf. 0), 

then the projection onto the last m factors 

pr : (x (0) ,x (1) , . . . ,x (m) ) G M n(m+1) i-> (x (1) , . . . ,i (n) ) G R nrn 

induces diffeomorphisms on the leaves of T . 

This can also be viewed as the fact that the foliation T corresponds to a connection 
A in the bundle pr : l n x ■•• x 1" = _> R« x • • • x W 1 = R nm with trivial 

curvature. The restriction of the projection pr to a leaf gives a one-to-one map. In 
this way there is a linear map among vector fields in ~R nm and (horizontal) vector fields 
tangent to a leaf. 

Note that the knowledge of this connection (foliation) gives us the superposition 
principle without referring to the function ^ (which we can change by composing 
it, for instance, with a diffeomorphism of M n ): if we fix the point X(o)(0), i.e. we 
choose a k = (k±, . . . , k n ), and m solutions x^(t), . . . , X( m )(t), then x^(t) is the unique 
point in M" such that (x( )(i), X(x)(t), . . . , X( m )(t)) belongs to the same leaf of T as 
(x( )(0), X(i)(0), . . . ,X( m )(0)). Thus, it is only T that really matters when the superpo- 
sition rule is concerned. 

On the other hand, if we have a connection V in the bundle pr : M n ( m+1 ) — > R nm 
with a trivial curvature, i.e. we have a horizontal distribution V in TM ra ( m+1 ) that is 



Superposition rules and Lie theorem 



6 



involutive, which therefore can be integrated to give a foliation in M"( m+1 ) such that 
the vector fields Y(t) belong to V (equivalently, are tangent to J 7 , i.e. are horizontal), 
then the procedure described above determines a superposition rule for the system (0). 

Indeed, let k G lR n enumerate smoothly the leaves JF fc of T (e.g. by a choice of a small 
cross-section of J 7 ), then . . . , £( m )(t); k) defined as the unique point X(o){t) G 

W 1 such that (x/q) (t) , x^ (t) , . . . ,X( m )(t)) G J-fc, is a superposition rule for the system 
(P) of ordinary differential equations. To see this, let us observe that the inverse is 
*(x (0 )(t),a;(i)(t), . . . ,X( m )(t)) = k, which is equivalent to (x {0) (t), x {1) (t), . . . , x {m) (t)) G 
Tk- If we fix k and take solutions x^(t), . . . , X( m )(t) of ((l}, then £(o)(£) defined by 
the condition \I/(x(o)(t), £(i)(t), . . . ,X( m )(t)) = k satisfies ((TJ). Indeed, let x', Jt) be the 
solution of (JTJ) with initial value 2:^(0) = X( )(0). Since the t-dependent vector fields 

Y(t) are tangent to J 7 , the curve (x', Jt),X(i)(t), . . . ,X( m )(t)) lies entirely on a leaf of 
J 7 , so on jF fc . But the point of one leaf is entirely determined by its projection pr, so 
£( )(t) = X(p)(t) and aJ(o)(t) is a solution. Thus we have proved the following geometric 
characterization of superposition rules: 

Proposition 1 Giving a superposition rule (p|) for a system of differential equations 
(PP equivalent to giving a zero curvature connection in the bundle pr : ]R( m+1 ) n — > R nm 
/or which the diagonal prolongations Y(t) of the t-dependent vector fields Y(t), t G M, 
defining the system (QJ) are horizontal. 

Note that we can also consider arbitrary manifolds iV instead of M n . The superposition 
functions are then given by maps $ : N m+1 — > N or by appropriate foliations in iV m+1 , 
i.e. zero curvature connections in the bundle pr : N m+1 — > N m . 

Remark 1 Frankly speaking, the connection is defined only generically, usually over 
an open-dense subset. But this is a general problem with superpositions, which hold 
only 'generically'. In the sequel all objects and constructions will be 'generic' in this 
sense. 

Example 1 Consider the (generalized) foliation T of codimension one generated by 
the vector field xd/dx + y d/dy on R 2 -{(0,0)}. It defines a zero curvature connection 
for the bundle pr : Rx(M-{0}) -> M-{0}, pr(x,y) = y el The leaves of J 7 are of the 
form (e t x,e t y), t G R, y ^ 0. In particular, the function ty(x,y) = x/y is constant on 
the leaves. The diagonal prolongation of the t-dependent vector fields a(t) x d/ dx are of 
the form a(t) {xd/dx + y d/dy) and are tangent to J 7 . This gives us the superposition 
rule for the linear differential equation x = a(t) x as follows: If x/u(t) is a solution, then 
(x(o)(£), X(i)(t)) belongs to one leaf, e.g with x/y = k, if and only if X(o)(i) = 
This is the standard superposition rule for this equation: 

$(x {1) (t);k) = kx (1) (t). 

Remark 2 It is clear from the above proposition that when the diagonal prolongations 
Y(t) generate a foliation of dimension smaller than the dimension of J 7 , we can change 
T (i.e. we can change the connection V) respecting the required properties. This means 
that the corresponding system admits many different superposition rules even if we 
regard compositions of $ with diffeomorphisms of M n as equivalence relations. 
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3 Lie theorem for ODE's systems admitting super- 
position rules 

In this section we shall give a proof of the classical Lie's theorem on ODE's admitting 
a superposition rule (cf. [TBI 121 131 EE]). 

First of all, in this proof we shall fill in all gaps which appear in the known lit- 
erature. Second, we shall use the alternative characterization of superposition rules 
given in Proposition 1 and last, but not least, the proof, as it will be presented, can be 
immediately extended to the case of partial differential equations. 

Remark first that all considerations will be 'generic' and local. Let us recall that a 
foliation T of an open dense subset of N = N m+l = N x ■ ■ ■ x N ((m+1) factors) defines 
a superposition rule for the system of differential equations (0), if and only if T is of 
codimension n, the projection pr:#->iVx---xA r = N rn (only m factors) onto the 
last m arguments maps leaves of T diffeomorphically and, furthermore, the generalized 
foliation T§ generated by the family {Y(t) \ t £ M} of diagonal prolongations of Y(t) is 
contained in T . 

We shall work only with the regular part of Tq. Such regular part is spanned by 
{Y(t) | t £ M.}, i.e. in any case, it is spanned by diagonal prolongations of some 

vector fields on N as [Y(t),Y(t')} = [Y(t),Y(t% etc. Let X u ...,X r be diagonal 
prolongations spanning locally the regular part of of dimension r, and therefore 
pr, (Xi ),..., pr^(X r ) are assumed to be linearly independent at a generic point. We 
clearly have r < mn. 

Since Xi, . . . ,X r span an r-dimensional foliation, then 

r 

[X a , Xp] = ^ 7 , 

7=1 

for some r 3 functions c a /j 7 defined on N. Note also that [X a , Xp] are diagonal prolon- 
gations as brackets of diagonal prolongations and the projections pr^X^, . . . , pr^pC,.), 
are assumed to be functionally independent. Then we shall use the following lemma. 

m 

Lemma 1 Let X a = ^^X Q ( a ) ; for a = 1, . . .r, and with r < mn, be diagonal pro- 

a=0 

longations to N of vector fields X a on N, and such that at each point p of N m the 

m 

vectors that are the projections of their values, pr jf (X a )(p) = '^^X a ^(p) , are linearly 

a=l 

r 

independent. Then, ~^^b a X a , with b a £ C°°(N), is again a diagonal prolongation if 

a=l 

and only if the coefficients b a are constant. 
Proof.- Let us write in local coordinates 

i=l 
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which implies that 

n m q 



Then, 



dx \a) 



n m r\ 

a 



^2b a (x {0) ,...,x {m) )X a = ^2^2^2b a (x i0) ,...,x {jn) )A % a (x ia) ) ^ , 

q=1 o=l i=l a=0 ( a ) 



which is a diagonal prolongation if and only if there are functions B l a (x), for a = 
0, . . . , m, and i = 1, . . . , n, such that for each pair of indexes i and a, 

r 

J^6q(x( ),. • .,X( m) )A l a (x {a) ) = B l a (x {a) ) , a = 0, ...,m % = l,...,n. 

a=l 

In particular, the r functions & a (x(o), . . . , £( m )) solve the following subsystem of linear 
equations in the unknown u a , for a — 1, . . . , r: 

^u a v4^0 (a) ) = £*(x( a )) , with i = 1, . . . ,n, a = 1, . . .,m. (11) 

a=l 

But the rank of the matrix (A l a (x( a )))a a is r < mn as the projections pr^Xi), . . . , pr !f (X r ) 
are assumed to be linearly independent. Thus the solutions u±, . . . , u r of (fTT| are unique 
and are completely determined by the matrix A^x^) 1 ^ and the vector B l a (x( a )) ha , with 
a — 1, . . . , m, so they do not depend on X( ). But since the diagonal prolongations are 
invariant with respect to the symmetry group S m+ i acting on N = N m+1 in an ob- 
vious way, the functions b a (x^, . . . ,x^) do not depend also on the other variables 

■ ■ ■ ,X( m ). 

Remark 3 Let us note that the assumption on the projections is crucial for the above 
lemma and actually without such assumption the result of this fundamental lemma is 
simply wrong as the following example shows. Consider the following two vector fields 
in R 2 which are prolongations of vector fields in IR: 

~ d d ~ d d 

ox ay ox oy 

and the functions 

bx(x,y)=xy, b 2 (x,y) = -(x + y) 

for which 

bi{x, y)Xi(x, y) + b 2 {x, y) X 2 (x, y) = - [x 2 — + y 2 — 



dx dy / 

is also a prolongation. However the coefficients 61 and b 2 are not constant. This is the 
standard gap in the proofs of Lie's theorem we found in the literature. One usually 
claims that a functional combination of diagonal prolongations is a diagonal prolon- 
gation only if the coefficients are constant without assuming that the corresponding 
projections are linearly independent. 
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Now, we can now prove the above mentioned Lie theorem using the previous results. 

Theorem 1 The system (QJ) on a differentiable manifold N admits a superposition rule 
if and only if the t-dependent vector field Y(t, x) can be locally written in the form 

r 

Y(t,x) = J2b a (t)X a (x) 

a=l 

where the t-dependent vector fields X a , a — 1, ... ,r, close on a finite-dimensional real 
Lie algebra, i.e. there exist r 3 real numbers c^ 7 such that 

r 

[X a ,Xp\ = ^c a/3 7 X 7 , Va,/3 = 1, . . . ,r. 

7=1 

Proof.- Suppose that the system admits a superposition rule and let T be the folia- 
tion corresponding to the superposition function. We know already that the generators 
{X a | a = 1, . . . , r} of the regular part of JF c T close on a Lie algebra 

r 

[X a ,X f3 ] = J2^f3" , X 1 , (12) 

7=1 

where the coefficients c Q/ g 7 are constant, so also 

r 

[X a , Xp\ = c a p 7 X 7 . 

7=1 

Since every Y{t) is tangent to Tq there are functions b"(x(o), ■ ■ ■ ,X(m)) such that 

r 

Y(t) = J2b?X a . 

a=l 

But Y(t) is a diagonal prolongation, so, using the fundamental lemma once more, we 
get that the b^ = b a (t) are independent on X(p), ■ ■ ■ , X{ m ). Hence 

r 

Y(t) = J2b a (t)X a (13) 

and also 

r 

Y(t) = Y,b a (t)X a . (14) 

a=l 

To prove the converse property, assume that the t-dependent vector field Y (t, x) can 
be written as in (|T3*j) and define Y{t) by (fTlj) . We can additionally assume that the 
vector fields X a are linearly independent over M. Thus they define an r-dimensional 
Lie algebra with structure constants c a/ 3 7 (and the corresponding simply connected Lie 
group if the vector fields are complete). 
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Since only non-trivial functional dependence of X\, . . . , X r is possible, there is a 
number m < r such that their diagonal prolongations to N m = N x • • • x N (m times) 
are generically linearly independent at each point. The distribution spanned by the 
diagonal prolongations X\, . . . , X r to N — N m+l is clearly involutive, so it defines an 
r-dimensional foliation JF of N. Moreover, the leaves of this foliation project onto the 
product of the last m factors diffeomorphically and they are at least n-codimensional. 
Now, it is obvious that we can extend this foliation to an n-codimensional foliation 
T with the latter property, and this foliation, according to proposition 1, defines a 
superposition rule. Here, of course, extending of a foliation means that the leaves of 
the smaller are submanifolds of the leaves of the extension. 

Remark 4 There is another way to look at the superposition the way we proposed. 
We can consider two projections: the first one, pr : N = N n ^ m+1 ^ — > fqnm j g on ^ e last 
m factors and the second one, pr : : N — > N the projection on the first factor. These 
projections are clearly transversal and the first one is a fibration if the foliation T is 
conserved, i.e. every curve in N m has a unique lifting in a fixed leaf of T . Then, the 
superposition associated with this lift is just pr x of the lift of particular solutions. 

Remark 5 We hope that it is clear to the reader that in our picture of superposition 
rules there is no real need to take all the manifolds N equal in the product where the 
superposition foliation lives. We can consider as well N — N x • • • N m with analogous 
projection and fibration property. This means that we can get a solution of a system 
on No out of solutions of some systems on N a , a = 1, . . . , m. This is a geometric picture 
for the Darboux (Backhand) transformations. We will, however, discuss this problem 
in a separate paper. 

4 Determination of the number m of solutions of a 
fundamental set 

Our proof of the Lie's Theorem contains an information about the number m of solutions 
involved in the superposition rule. For a Lie system defined by a t-dependent vector 
field of the form 

r 

Y(t,x) = ^2b a {t) X a {x) , 

a=l 

with generic b a (t) the number m turned out to be the minimal k such that the diagonal 
prolongations of X\, . . . , X r to N are linearly independent at (generically) each point: 
the only real numbers solution of the linear system 

r 

^c a X Q (x (a) ) = 0, a = l,...,fc 

a=l 

at a generic point (xm, . . . , £(&)) is the trivial solution c a — 0, a — 1, . . . , m, for k — m 
and there are nontrivial solutions for k < m. 
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For instance, for the Riccati equation 

x = b (t) + b±(t) x + b 2 {t)x 2 , 

in which the vector fields generating the foliation JF are prolongations of the vector 
fields 

Xo = ^~, Xi = x-?~, X 2 = x 2 ^-, (15) 
ox ox ox 

representing an sl(2, R)-action, we see that the system 

C + C\X\ + c 2 x\ = , c + c x x 2 + c 2 x\ = 
has a nontrivial solution but the one given by 

C + C\X\ + C 2 x\ = , C + C X X 2 + C 2 x\ = , C + C1X3 + c 2 x\ = , 

does not admit non-trivial solutions because the determinant of the coefficients is in- 
vertible when the three points x±, x 2 and £3 are different. This implies that m = 3 in 
the superposition rule for the Riccati equation. 



5 Nonuniqueness of the superposition rule 

In some cases the foliation JF spanned by the prolongations of t-dependent vector fields 
defining the dynamics is already n-codimensional and we get a unique (minimal) su- 
perposition rule (e.g. Riccati equation). In the cases with codimjF > n, we have some 
ambiguity in choosing the superposition rule as we can extend JF to an ra-codimensional 
foliation in different ways. 

Example 2 Consider the action of the Abelian group K acting on N = M 2 by horizontal 
translations., i.e. X\ = d/dx. This action on IR 2 is free and we have m — 1, so 
N = 1R 2 x M 2 and To is spanned by d/dx^ + d/dx/u, where the coordinates in N are 
denoted (x( ), 2/(0), £(1), V(i))- We can extend such foliation to a 2-dimensional foliation T 
of N with the required property with respect to the projection pr(x( ), j/(o), j/m) = 
(x(i), in different ways. For instance, we can take JF to be given by the level sets 
of the mapping 

F{x(o), 2/(0), &(i), £/(i)) = (x (0 ) - /(2/( ), ?/(i))) 

with / being an arbitrary function such that df/dy^ 7^ 0. Then, every solution 
(x(i)(t), U(i){t) = 2/(1) (0)) of the system of differential equations 

x = a (t) , y = , 

gives a new solution (x(o)(t), 2/(o)(^) = 2/(o)(0)) associated with the level set of (fci, /c 2 ) by 



(x (0 )(t) = x ( i)(t) + fci, 2/(0) (t) = 2/( )(0)) 
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where 2/(o) (0) is the unique point in R satisfying 

f(y {0) (0),y {1) (0)) = k 2 . 
In the case f{y(o),y(i)) = V(o) — V(i) we recover the 'standard' superposition rule: 

3/(1)1 fo) = ( x (i) + h, + fa) ■ 

Example 3 A very simple example is given by the separable first-order differential 
equation 

x = a (t) f(x) , 

with a and / being arbitrary smooth functions, and where / is assumed to be of a 
constant sign (otherwise we can restrict ourselves to a neighbourhood of a point in 
which / does not vanish) . In this case N = R and we can consider the vector field in R 

*(*) = /(»)|. 

As the function / does not vanish, we have m = 1 and the diagonal prolongation 

d d 

X(Z( ),Z(1)) = /(X( 0) ) 77— + Ml)) 77— 

OX(o) OX(i) 

generates a one- dimensional foliation in R 2 whose leaves are the level sets of a function 
\P(x(o), such that 

/( x (o)) 7^ + /0(i)) 7^ = , 

dx(p) dx {1) 

which gives rise to the following characteristic system 

dx(o) dx^i) 

Therefore, if the function 4>(y) is defined by 

m = [jk)' 

then we find that the leaves are characterized by a constant k in such a way that 

0(«(o)) - <f>(x(i)) =k. 

The function is a monotone function, because <fr'(x) = f(x) and f(x) has constant 
sign. Therefore, there exists an inverse function which allows to write the superposition 
rule as 

x = 0" 1 (k + 0(X(i))) . 

For instance, if f(x) = 1/x 2 , we find that <f)(x) = —1/x = an we obtain the 

following superposition rule. 

x - 



1 — kx(\) 
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Example 4 It has been pointed out in the following example of the original Lie's 
work: 

— = a 12 {t)y + b 1 {t) 
ay 

— = -a 12 (t) x + b 2 (t) 

In principle, it is a particular example of an inhomogeneous linear system and we expect 
to have an affine superposition rule involving three different solutions: 

x = $i(x(i),a; ( 2),X(3)) = x {1) + k x (x {2 ) - + ^0(3) - £(i)) • 

However, we can obtain a superposition rule which is not linear but involves only two 
solutions. It is due to the fact that here is not the affine group in two dimensions which 
is play the relevant role, but the Euclidean group. The foliation corresponding to this 
Lie system is generated by the prolongations of the vector fields 

v 9 v 9 v 9 9 

ox Oy ox oy 

First of all, m is different from 1, because there exist nontrivial coefficients Ai, A2 and A3 
such that Ai X x {x^) + \ 2 X 2 {x^) + A 3 X 3 (x^) = 0, at a given point x^, for instance, 
Ai = —y(i), Ai = A 3 = 1. However, the only coefficients Ai, A 2 and A 3 such that 

\iXi(x( 1) )+\ 2 X 2 (x( 1) )+\ 3 X 3 (x(i)) = 0, AiXi(x(2))+A2X2(x(2))+A 3 X3(x( 2 )) = 0, 

with 7^ X( 2 ) are Ai = A 2 = A 3 = and therefore m = 2. 

The function \1/ defining the superposition rule satisfies Xi^ = X 2 ty = X 3 ty = 
with 

~ d d d ~ d d d 

Xl = Ti — ^ Ti ^ Ti — ' X 2 — — ~\~ — h — — , 

ox 0x1 ox 2 Oy oyi oy 2 

and 

~ d d d d d d 

X 3 = y-£- - x— + yi- xi— + y 2 - x 2 — . 

ox Oy 0x1 Oyi ox 2 OyZ 

The two first conditions imply that \& must be of the form 

*(sc , 2/o, ari, yi, x 2 , 2/2) = *P(xo -xi,x - x 2 , y - y u y - y 2 ) , 
what suggests the change of variables 

ui=x -xi, u 2 = x -x 2 , u 3 = x 0l i>i = y -yi, v 2 = y - y 2 , v 3 = y , 
and then the third condition X 3 ty = is written 

d d d d 
v ±x ^ v ^a «i q M 2^— , 

OUi OU 2 OVi ov 2 

for which the characteristic system is 

dui du 2 dvi dv 2 



Vi V 2 -Ui —u 2 
from where we find the first integrals 

ul + vf = (x - xi) 2 + (y - yi) 2 =C X , u\ + v\ = (x - x 2 f + (y„ - y 2 ) 2 = C 2 

which determine the superposition foliation and provide us with a superposition rule 
for the given system involving only two particular solutions (i.e. with m = 2). 
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6 Lie systems in Lie groups and homogeneous spaces 

Let us consider now the particular case m — 1, i.e. when a single solution is enough 
to obtain any other solution. Let us assume additionally that T = J-"q, i.e. that the 
superposition rule is unique as a foliation. This means that r = n and the vector fields 
Xi, . . . ,X n generically span QTN. Assume for simplicity that they span QTN globally 
and are complete vector fields. Since these vector fields close on an n-dimensional Lie 
algebra q, we conclude that there is a transitive action on N of the simple-connected n- 
dimensional Lie group G associated with q, so that N = G/H with a discrete subgroup 
H of G and the foliation JF = JF is generated by the fundamental vector fields of the 
G-action. If H is trivial and we consider the standard action of G on itself by left 
translations L g , the vector fields Xi are just right-invariant vector fields on G. As a 
superposition function corresponding to T = T§ we can choose the group multiplication 
$ : G x G — > G, $(g(i), k) = 9(i)k- m this case \I/(g(o), g(i)) = g7^ is left invariant 

y(g'9(o),g'g(i)) = *(#(o),#(i))- 

Conversely, given a Lie system defined by a t-dependent vector field of the form 

r 

Y(t,g) = J2Ut)X R a (g), 

a=l 

where X* is a basis of right-invariant vector fields, then the projectability condition is 
satisfied and there is a uniquely defined superposition rule. Note however that if the 
vector fields X* generate a smaller Lie subalgebra, the superposition rule is not unique. 

Let {ai, . . . ,a n } be a basis in T e G. This linear space can be identified with the 
Lie algebra g of G, the set of left invariant vector fields on G: for each a e T e G let 
X^ denote the corresponding left-invariant vector field in G, given by X^(g) = L gJfe a. 
Similarly, X* denotes the right-invariant vector field in G given by X*(g) = R g * e a - A 
curve in T e G, 

n 

a(t) = ^ b a (t) a a , 

a=l 

gives rise to a t-dependent vector field on G 

n 

X R (t,g) = X« a{t) (g) = Y J b a (t)X«(g) , (16) 

where in the right hand side X^ is a shorthand notation for X* a - The associated system 
of differential equations determining the integral curves of such a t-dependent vector 
field reads 

n 

and applying Rg-i( x )*g(x) to both sides of (fTTj) we find the equation 



n 

Rg-Ht^gi^git) = ^fra(£)a a = a(t) , 

a=l 



(18) 
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that with some abuse of notation we will write 



(jr)(t) 



a(t) . 



(19) 



The solution of this equation starting from the neutral element can be solved by making 
use of a generalization of the method developed by Wei and Norman ^SJ El] (see e.g. 
jH] and [3]) for solving an analogous linear problem. A simpler situation is when a(t) 
takes values not in the full Lie algebra g = T e G but in a subalgebra. 

Now, let H be an arbitrary closed subgroup of G and consider the homogeneous 
space N = G/H. Then, G can be seen as a principal bundle r : G — > G/H. Moreover, 
it is also known that the right-invariant vector fields are r-projectable and the r- 
related vector fields in N are the fundamental vector fields — X a = —X aa corresponding 
to the natural left action of G on N, T* g X*(g) = —X a (gH). In this way we can associate 
with the Lie system on the group X given by (fTH|) a Lie system on N: 



Therefore, the integral curve of (|20|) starting from x are given by x(t) = &(g(t),x ), 
with g(t) being the solution of (fTSj) with g(Q) = e. 

Let us note that even if the original t-dependent vector field on N is projectable to 
Nq by means of a submersion ir : N — > No, the superposition rule cannot be projected in 
general and the number m of solutions appearing in the superposition rule changes. For 
instance, in the case of Riccati equation a fundamental set is made of three solutions, 
while for the linear realization of SL(2, IR) on IR 2 only two solutions are needed, and for 
SL(2,M.) acting on itself only one solution is sufficient. 

In the particular case of a Lie system g g~ l = a in a Lie group G, it was recently 
shown [1] that the knowledge of a particular solution of the corresponding system in a 
homogeneous space for G reduces the problem to one on the isotopy group of a point in 
the homogeneous space. So, if x(t) is the particular solution of the associated Lie system 
in a homogeneous space starting from xq, then we can choose a curve g(t) such that 
<&(g(t),xo) = x(t) and there should exist a curve h(t) € G Xo such that g(t) = g(t) h(t). 
Such curve h(t) is a solution of the Lie equation in G Xo , hh" 1 = Ad(?~ 1 (a + gg^ 1 )- 
Therefore, finding a solution of such equation in the subgroup G xo , we can recover the 
solution g(t) of the Lie system in G as g(t) = g(t) h(t). Using a new solution starting 
from a new point, the problem is further reduced and therefore with a number of known 
solutions we can directly write the general solution. 

Another relevant case is when there exists an equivariant map F : N\ — > N2 between 
two homogeneous spaces of a Lie group G. In this case the corresponding fundamental 
fields are F-related and then the image under F of an integral curve of a Lie system in 
Mi is an integral curve of the corresponding system in M 2 . A very simple example is 
the following: the function F : R 2 - {(0, 0)} -> 1 = R U {00} given by 



71 




(20) 



a=l 




x 2 = 
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is equivariant with respect to the linear action of the Lie group SL(2, R) on M 2 — {(0, 0)}, 
and its action on the completed real line M, given by 

t / * \ ax + (3 .5 

${A,x) = — , ifx^ — , 

7 x + o 7 

oo) = -, -5/i) = oo, 

7 



when A is the matrix given by 



A= [ a P ) e SX(2,R) . (21) 



7 5 

Consequently, given an integral curve of the system of differential equations 

±i = ^b 2 x 1 + b 1 x 2 , 

x 2 = -b 3 xi--b 2 x 2 , 

the curve x(t) = xi(t) / x 2 (t) is an integral curve of the Riccati equation 

x = b\ + b 2 x + b 3 x 2 . 
This is precisely the method by which Riccati arrived to this last equation. 

7 Partial superposition rules 

Finding new solutions from know ones is a very usual method in the theory of both 
ordinary and partial differential equations and this procedure has many applications in 
physics. In order to deal in a geometric way with this problem we introduce next a 
concept generalizing those of nonlinear superposition introduced so far and the one in 
|T2"] as well as that of connecting function used in ^Hj (see also [Hj for some examples). 

A partial superposition rule of rank s of m solutions for the system of ordinary 
differential equations (0) is given by a function $ : ]R nm+s — > M. n , 

x = $(x ( i), . . . , x (m) ; fci, . . . , k s ) , (22) 

such that if {xr a -\(t) \ a — 1, ... , m} is a set of m particular solutions of the system (|Tj). 
then, at least for sufficiently small t, 

x(t) = $(x ( i)(t), . . . ,X( m )(£); fa, ...,k s ) , (23) 

is also a solution of the system (0), where k = (fa,...,k s ) is a set of s arbitrary 
constants. 

Note that this new concept reduces to the previously considered superposition rule 
for s = n and coincides with a s-parameter family of connecting functions in the sense 

of mg. 
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Given such a partial superposition function, there is a non-uniquely defined function 
* : BH m+1 ) -> M n such that 

T i / \ f fcj if i < s 

* (x (0) ,x (1) ,...,x (m) ) = | Q , f . > s 

The last (n — s) equations are restrictions defining a submanifold M of co dimension 
(n — s) of R n ( m+1 ) ; i.e. of dimension mn + s, and the other equations define a foliation of 
codimension s in M. Now, following the same procedure as in Section 2 we will arrive 
to a distribution in M n ( m+1 ) spanned by the t-dependent vector fields Y a defined by 
(110)1 which provide us with a distribution in M, because they are tangent to M, whose 
integral leaves are n m-dimensional and each leaf is fixed by the choice of s constants 
k±, . . . , k s . 

Moreover, the restriction pr| M of pr on the submanifold M defines a subbundle of 
pr : M n ( m+1 ) — > M. nm and establishes diffeomorphisms among the different leaves and 
allows us to identify among them the leaves of the foliation defined by the prolongation 
of the given non-autonomous system. 

Conversely, if M is a submanifold of R n ( m+1 ) of codimension (n — s) which defines a 
subbundle of pr : R n ( m +V -> R nm such that the distribution defined in M n ( m+1 ) by the 
prolongation of vector fields is also a distribution in M, i.e. such vectors are tangent 
to M, and the restriction pr| M provides us with diffeomorphisms among the different 
leaves allowing us to identify among them the leaves of the foliation defined by the 
prolongation of the given non-autonomous system. Such diffeomorphisms can be used 
to define a superposition rule of m solutions involving s constants, 

If for instance we consider the linear system 

^ = a n (t) x 1 + a 12 (t) x 2 (24) 
dx 2 

— = a 21 (t) x 1 + a 22 {t) x 2 (25) 

then it admits a superposition function of rank one and involving one particular solution, 
F(xny, k) = k im, which determines the three-dimensional subbundle M of pr : IR 4 — > 
M 2 defined by the restriction to the subset given by the relation 

(1) (1) ' 

which is endowed with a foliation: each leaf is characterized by a real number k and is 
defined on the set of points (x 1 , x 2 , zL, x 2 ^) such that x l x 2 ^ — x 2 a?jL = 0. 

However we have also a superposition function of rank one but involving two con- 
stants: 

F(x (1 ), 2(2); k) = 2(i) + k 2(2) , 

The subbundle now will be defined by 

X (2)( x2 — X (l)) ~ X f2)( xl ~ X (l)) — • 
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8 Superposition rules for PDE's 

Consider now the system of first-order PDE's of the form: 

— = Y:(t,x), xeR n , t = (t\...,t s )eW, (26) 

whose solutions are maps x(t) : IR S — * M. n . 

A particular case of (|26|) when s — 1 is |I} . The main difference of (}26|) with respect 
to is that for s > 1 we have no, in general, existence of a solution with a given initial 
value x(0) G M n . For a better understanding of this problem, let us put (f2T)|) in a more 
general and geometric framework. 

For a manifold N of dimension n consider the trivial fibre bundle 

P^ = R s x N -> W . 

A connection Y in this bundle is a horizontal distribution in TP^. i.e. an s-dimensional 
distribution transversal to the fibres. It is determined by horizontal lifts of the coordi- 
nate vector fields d/dt a in W which read 

Y a = ^- a +Y a (t,x) 

with 

Thus, the solutions of (j2B| can be identified with integral submanifolds of the distribu- 
tion Y, 

(t,Y(t)), teR s . 

It is now clear that there is an (obviously unique) solution of (I26j) for every initial data 
if and only if the distribution Y is integrable, i.e. the connection has a trivial curvature. 
This means that 



[Y a ,Y b ]=Y J fab C Y c 



c=l 



for some functions f a t, c in P^. But the commutators [F a ,y"h] are clearly vertical while 
Y c are linearly independent horizontal vector fields, so f a b c = which yields the in- 
tegrability condition in the form of the system of equations [Y a ,F fe ] = 0, i.e., in local 
coordinates, 

-£{t,x) - -^(t,x) + [Yi(t,x)-^(t,x) - Yi{t,x)^{t,x)\ =0. (27) 

j=l v / 

Let us assume now that we work with a system of first-order PDE's of the form (|26J) and 
satisfying the integrability conditions (|27j). Then we are sure that, for a given initial 
value, there is a unique solution of (j2f)jl . Now, we can think about superposition rules 
for such solutions. It is, however, completely obvious that the concepts of superposition 
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rules we have developed can be applied with no real changes to our case of PDE's. In 
the formula (J2J) we should now think that t is not a real parameter but i 6 I s . The 
only difference when passing to the foliation induced by the superposition function \1/ 
is that we differentiate (jHJ) not with respect to the simple parameter t but with respect 
to all t a . Therefore, the proposition 1 takes the form: 



Proposition 2 Giving a superposition rule for the system \2b\) satisfying the integrabil- 
ity condition \27\) is equivalent with giving a connection in the bundle pr : A r * m+1 ) — > N m 
with a zero curvature and for which the diagonal prolongations Y a (t) of all the vector 
fields Y a (t), t e W , a = 1, . . . , s, are horizontal. 

Also the proof of Lie theorem remains unchanged. Therefore we get the following analog 
of the Lie theorem for PDE's: 



Theorem 2 The system \2b}) of PDE's defined on a manifold N and satisfying the 
integrability condition \2T\) admits a superposition rule if and only if the vector fields 
Y a (t,x) on N depending on the parameter t e W , can be written locally in the form 

r 

Y a (t,x) = Y,<(t)XM, a = l,...s, (28) 

a=l 

where the vector fields X a close on a finite- dimensional real Lie algebra, i.e. there exist 

a/3 



r 3 real constants c n n 7 such that 



[X a , Xp\ — c a p 7 X^ . 

7=1 

Let us observe that the integrability condition for Y a (t, x) of the form (|2"K|) can be 
written as 

r 

q,/3,7=1 



X 7 = 0. 



Example 5 Consider the following system of partial differential equations on M 2 asso- 
ciated with the sl(2, M)-action on M represented by vector fields f)15|) : 

u x = a(x,y)u 2 + b(x,y)u + c(x,y) , 
u y = d(x,y)u 2 + e(x,y)u + f(x,y) . 

This equation can be written in the form of a 'total differential equation' 

(a(x, y)u 2 + b(x, y)u + c(x, y))dx + (d(x, y)u 2 + e(x, y)u + f(x, y))dy = du . 

The integrability condition just says that the one-form 

to = (a(x, y)u 2 + b(x, y)u + c(x, y))dx + (d(x, y)u 2 + e(x, y)u + f(x, y))dy 



is closed for arbitrary function u = u(x,y). If this is the case, then there is a unique 
solution with the initial condition u(xq, yo) = uq and there is a superposition rule giving 
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a general solution as a function of three independent solutions exactly as in the case of 
Riccati equation: 

(U - M(l))(M(2) - U(3)) = , 
(«-U( 2 ))(«(i)-U(3)) 

or 

^ _ ('»(!) ~ U {Z) )u [2) k + M(i)(tt( 3 ) ~ M( 2 )) 
(«(1) - «(3))fe + («(3) - "(2)) 

9 Concluding remarks 

In this paper we have identified and solved a gap present in previous proofs of the 
necessary and sufficient conditions for the existence of a superposition rule for ordinary 
differential equations. In doing this we provided the superposition rule with a much 
better geometrical interpretation which allows us to exhibit many interesting properties. 
For instance, it is now clear by inspection that the vector field as given in the Theorem 
1 may be multiplied by a function of time, i.e. we may perform a re-parametrization in 
time, and still get an equation admitting a superposition principle. In this way we find 
that if the superposition rule holds true for a vector field field, it is also true for all re- 
parametrized vector fields. In particular, it would allow us to write a superposition rule 
also for vector fields which may be reduced to autonomous ones via re-parametrization. 
The new geometrical interpretation also paves the way to a proper treatment of the 
superposition rule for partial differential equations. We hope to be able to extend the 
treatment to field theory and perhaps be able to get interacting field theories out of free 
ones very much as it happens for ordinary differential equations. Indeed it is known that 
Riccati equation obtains from a linear system. In previous papers it has been shown how 
to cast completely integrable systems in a generalized version of Lie-Scheffers systems, 
as most completely integrable systems do arise as reduction of simple systems we hope 
to be able to show in general that systems which allow for a superposition rule may be 
derived from 'simple ones', both for ordinary and partial differential equations. 
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